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The indirect vector control system can realize the on-line identification of 
the rotor time constant based on electromagnetic torque. However, still 
effective range of the torque model is suitable only for heavy load 
conditions. In this paper, an improved torque tracking model is constructed 
by modifying the conventional torque model with rotor time constant 
identification. Through the establishment of the small signal model of the 
identification system, the stability of the torque model is deeply analyzed, 
and it is concluded that the torque model has two balance points, namely the 
expected balance point and the error balance point; the error balance point 
limits the scope of the torque model root cause. The simulation show that the 
problem of unbalance point can be avoided with the aid of the proposed 
algorithm. In addition, the design ideas of proportional-integral adaptive law 
parameters are given from the two perspectives of system stable operation 


Tracking and optimal convergence process. 
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1. INTRODUCTION 

The indirect field oriented control (IFOC) scheme is widely used in high-performance induction 
motor drive occasions due to its simple structure and high reliability [1]. The rotor time constant is a key 
parameter to achieve rotor field orientation and high-performance control [2]. However, in practical 
applications, the rotor inductance will change with the degree of magnetic saturation of the motor, and the 
rotor resistance will be affected by changes in the internal temperature of the motor [3], [4]. The deviation of 
the rotor time constant will lead to wrong field orientation, which will deteriorate the system control 
performance [5], [6]. 

Many researches have conducted in-depth research and proposed many online identification 
schemes for rotor time constants. Hasan and Husain [7] combines the advantages of Romberg observer and 
sliding mode observer, and proposes a Romberg-sliding mode observer to realize on-line identification of 
rotor resistance. Tran et al. [8] uses dual sliding mode surfaces to achieve simultaneous correction of rotor 
resistance and slip gain (the reciprocal of the rotor time constant). Ammar et al. [9] uses the voltage response 
information in the transition process to modify the rotor time constant. In [10], [11] uses the high frequency 
information of certain physical quantities in the reference adaptive system model (RASM) speed observer to 
realize the identification of the rotor time constant. In [12]-[15] uses signal injection, artificial neural 
network, and extended Kalman filter to realize the on-line identification of the rotor time constant. Among 
many identification algorithms, the identification algorithm based on RASM is favored by people because of 
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its easy implementation and good stability [1]. Many algorithms has different implementation forms 
depending on the selected physical quantity, such as q-axis rotor flux model [16], [17], d-axis stator voltage 
model [18], [19], reactive power model [1], [5], [20], [21], electromagnetic torque model [22]-[25]. 
Abdalla et al. [20] established a unified mathematical model of this type of identification system, revealing 
the internal connection between different identification models. Bermúdez et al. [23] comprehensively 
compared several identification schemes from the perspectives of deviation sensitivity, convergence speed, 
parameter dependence, and pointed out that each scheme has its limitations; according to different 
applications, the optimal identification model can be selected. Rashag et al. [21] pointed out through steady- 
state analysis that the rotor time constant identified by the torque model cannot converge to its true value 
under certain load conditions. Mohamed and Habbi [22] pointed out that in the full load interval, the polarity 
of the torque deviation is not always consistent. Bermúdez et al. [23] proves through experiments that when 
different rotor time constants are deviated, the torque-current curves overlap in the light load region. 
Cao et al. [24] uses the unified RASM theory to prove that the torque model does not satisfy the lyapunov 
stability condition when the load is small. The existing In has a more comprehensive analysis of the torque 
model, and pointed out that the torque model is only suitable for the heavy load condition of the motor. 
Nevertheless, the research on the torque model still has the following shortcomings: the stability of the 
identification system requires further in-depth analysis; for the problem of the limited action interval of the 
torque model, the existing in rarely mentions the improvement scheme. Chutisant [25], proposed interesting 
extended particle filter to recognition hand writing. Theparcticl filter can used as a new idea that may be 
directly applied to the stability analysis of the rotor time constant identification system based on the torque 
model. 

In this paper an improved torque model by simply modifying the torque tracking equation in the 
conventional torque model was proposed. There are two balance points in the torque model, namely, the 
desired balance point and the error balance point by solving the steady-state balance point. The value of the 
rotor time constant at the unbalance point not only depends on the real rotor time constant value of the motor, 
but also is affected by the load ratio (the ratio of the torque current to the excitation current). The model can 
avoid unbalance point is the leading factor that causes the instability and limitations of the torque model, and 
the selection of inappropriate regulator parameters will cause system instability. In addition, the model 
regulate parameters from two perspectives: the stable operation of the system and the optimization of the 
convergence process. 


2. THEORETICAL BCKGROUND 
2.1. Mathematical model of induction motor 

In the dq synchronous rotating coordinate system, the dynamic mathematical model of induction 
motor based on indirect vector control can be expressed as [2]: 


F dad. š L Lm d 
Vs = Reis + OLs is +Jz@c0Lsis + Jawe Tbr +T uYr (1) 
d L al 
avr T T Dr Tap” — hosp (2) 
— 5a 
Ws = (3) 


Bok 
Trloq 


where: Vs=[VsaVsq] > is=[isaisg]’. Waa sql” are the stator voltage and stator respectively in the dq 


synchronous rotating coordinate system Current and rotor flux linkage; t, is the real rotor time constant of 


Hon 


the motor, represents the estimated value; @e, @,p, @, are the synchronous electrical angular velocity, 


rotor electrical angular velocity, and slip electrical angular velocity, respectively; R, is the stator resistance, 
Lm» Ly. Ls are the motor excitation inductance, stator inductance, and rotor inductance respectively; Jz = 


h A oiim am superscript "*" indicates the command value. 
1 0 LsLr 
2.2. Conventional torque model 

The torque model proposed is realization block diagram [22]. The actual torque expression has been 
incorporated into the adaptive PI gain. The mathematical expression of the conventional torque model can be 
expressed as [22]: 
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Teref = ispWra = isaWrp = isqWra = isaWrq (4) 
Teest = Lmisqisa (5) 
=+ kp AT + f kirATdt (6) 


where: AT -Teref Teest) Kpr, and kir are proportional and integral adaptive gains respectively. Due to the 


symmetry of the torque expression with respect to the torque current, the following will take the positive 
torque current condition as an example for analysis. 


3. PROPOSED METHOD 
Figure 1 shows the block diagram of the rotor time constant online identification system based on 
the improved torque model. As shown in Figure 1, isq Y,a isd Y, a is selected as the physical quantity in the 


traditional torque model. In the case of accurate magnetic field orientation, that is, y, ,=Lmisq and y, q™® the 


physical quantity is equal to Lmisqisa: 


Xref = isqWra = isaWrq (7) 
Xest = Lmisqisa (8) 
A: aE 

A kpAX + f kirAXdt (9) 


where AX=X;e¢ Xesti Kp» Ki are proportional integral coefficients. Let the differential term in (2) be 0, the 
steady-state rotor flux in the dq coordinate system can be obtained as (10). 
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The steady-state AT expression in the conventional torque model can combine (4), (5) and (10) and be 
obtained (12). 
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It can be seen that compared with the conventional torque model comparing (11) and (12), the estimated slip 
gain value at the error balance point EP; of the proposed scheme is a negative number -1/(z,k;’). 
Therefore, the proposed algorithm can naturally avoid the problem of wrong balance point. If the stability of 
the system at the desired balance point is also ensured, the proposed scheme is globally stable. 

The RASM method construct with a new physical quantity X = i,gy,, t isd y, g In the case of 
accurate magnetic field orientation, the physical quantity mentioned is also equal to Lmisqtsa. This is why the 
algorithm proposed in this paper is called an improved torque model. In fact, the variable X is a fictitious 
physical quantity with no physical meaning. Using the variable X as the set physical quantity in the RASM 
can naturally avoid the unbalance point in the traditional torque model, which will be described detail in 
section result and discussion. 
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Figure 1. Overall proposed model 


4. RESULTS AND DISCUSSION 

In this section the conventional torque model with rotor time constant identification models were 
developed for 3hp induction motor. In order to verify the correctness of the proposed theory and design 
scheme, a simulation platform composed of a set of 10 kW induction motor as shown in Table 1. The table 
consists in the rated value of power, voltage, curremt, frequency, stator and rotor resistance, rotor nd stator 
inductance. 


Table 1. Induction motor parameters 


Parameter Value Parameter Value 
Power (kW) 10 Stator resistance (Q) 0.076 
Voltage(V) 220 Rotor resistance (Q) 0.055 
Current(A) 32 Stator inductance (mH) 14.1 
Frequency (Hz) 100 Rotor inductance (mH) 14.1 
Pairpole 2 Mutual inductance (mH) 13.6 


4.1. Conventional torque model experiment 

In this section also investigate the influence of PI regulator parameters on the stability of the 
identification system [18]. The experiment in this section will focus on verifying the correctness of the 
theoretical analysis. Figure 2 shows the identification simulation waveforms of the torque model with 
different PI parameters under the conditions of i,, = 40 A and n=1200 r/min. According to (7), it can be seen 
from Figure 2(a) that under the condition of kpr=0, the integral coefficient must satisfy kir <2.8 to ensure 
the stability of the system. From Figure 2(b), when kir = 4, it can be seen that the identification system 
collapses. According to (6), it can also be obtained: when the integral coefficient is fixed at 1, the 
proportional coefficient must be less than 0.721 to avoid instability. From Figure 2(c), it can be seen that 
when k,7=0.8 and kir = 1, the identification system collapses. It can be found that the simulation results are 
consistent with the theoretical analysis, which proves the correctness of the theoretical analysis. 
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Figure 2. Simulation result of rotor time constant identification of conventional torque model: (a) Kp7=0, 
kiT=1, (b) Kpr=0, kiT=4, and (c) Kpr=0.8, kiT=1 


4.2. Proposed torque model 

When the motor speed is 1200 r/min, this section selects two typical torque current values for 
experiments to verify the proposed algorithm, is;,=10 A and 40 A. When k;=1,3 Figures 3(a) and (b) 
respectively show the rotor time constant identification simulation waveforms of the proposed RASM under 
two working conditions. From Figure 3, it can be seen that the proposed algorithm can achieve accurate 
identification of the rotor time constant under different torque currents, which proves that the proposed 
algorithm can effectively broaden the working range of the conventional torque model. It can be found from 
the simulation results: when the torque current is small (10 A), with the increase of k;, the identification 
parameter convergence speed will increase; when the torque current is large (40 A), with the increase of k; 
Increase, the identification parameters appear under-damped oscillation process. Under the same two groups 
of torque currents, Figure 4 shows the results of the rotor time constant identification experiment using the 
proposed integral coefficient design scheme for case (a) at isg=10 A and for case (b) ig=40 A. From 
Figures 4(a) and (b), it can be found that the optimal parameter convergence process can be realized at 
different operating points using the design strategy proposed in this paper. 
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Figure 3. The rotor time constant with different integral coefficients: (a) isg=10 A and (b) ig=40 A 
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Figure 4. The online rotor time constant of the integrated coefficient: (a) isg=10 A and (b) isg=40 A 
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Figure 5 shows the rotor time constant identification experiment results when the motor runs at 
30 r/min and torque current 40 A. From Figure 5, it can be found that the proposed algorithm can also 
achieve accurate identification of the rotor time constant at lower speeds. When the motor speed is fixed at 
1200 r/min, the torque current command changes according to 40, 30, 20, 10, 40, and 10 A. At this time, the 
corresponding motor phase (A) current are shown in Figure 6(a), and the slip gain identification waveform is 
being presented in Figure 6(b). From Figure 6, it can be seen that the identification value of the rotor time 
constant changes with the torque current command. The reason for this phenomenon is the existence of non- 
ideal factors such as iron loss and stray loss in real motors. 
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Figure 5. The online identification of the time constant of the RASM rotor 
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Figure 6. The slip gain and i,, during torque current step: (a) i,,simulation result and (b) differential gain 


Table 2 shows the actual output torque measured by the torque meter after the deviation of the 
theoretical output torque, the rotor time constant for orientation, and the identification algorithm for different 
torque current commands. The simulation result is that the motor is tested at a stable speed of 1200 r/min. 
From Table 2, it can be found that when 1/7, =8s—1 and the torque current command value is small, the 
actual output torque is greater than the theoretical torque. As the torque current increases, the actual output 
torque will be less than that. Theoretical value. When the identification algorithm is enabled, the torque 
control accuracy can be effectively improved. Table 2 shows, from the simulation data it also can observe 
that even after the identification algorithm is enabled, the output torque is still slightly less than its theoretical 
value. The reason for this phenomenon is that the motor iron loss branch will occupy a part of the stator 
current g axis component. How to compensate the impact of iron loss and other non-ideal factors on the 
torque control accuracy will be the focus of the author's subsequent research work. 


Table 2. Theoretical and measured torque values before and after identification 
Tracking torque (N-m) 
l/tr=8 After fixing 


Current (A) Reference torque (N-m) 


5 3.79 6.32 3.23 
10 7.68 9.24 6.91 
15 11.6 10.61 10.82 
20 15.2 11.91 14.25 
25 19.3 12.92 18.32 
30 23.1 14.52 22.01 
35 26.8 15.61 25.52 
40 31.2 17.59 29.57 
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5. CONCLUSION 

This paper introduced an improved torque model by modifying the reference torque calculation 
equation. The model integrate two indicators of stable operation of the system and optimized convergence 
process to give the idea of selecting the PI parameters. The experiment verifies the correctness of the 
theoretical analysis and the proposed scheme. In order to further improve the motor control performance, 
future work will focus on the estimation of the rotor flux linkage and the compensation of non-ideal factors 
such as iron loss. 
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